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SUMMARY 


Continuous fiber silicon-carbide/aluminum composite laminates 
with slits were tested statically to failure. Five different layups 
were examined: [0] 8 , [0 2 /±45] s , [0/90] 2s , [0/±45/90] s , and [±45] 2s . 

Either a 9.5 or a 19 mm slit was machined in the center of each 
specimen. The strain distribution ahead of the slit tip was found 
experimentally with a series of strain gages bonded ahead of the slit 
tip. A three-dimensional finite element program (PAFAC) was used to 
predict the strain distribution ahead of the slit tip for several 
layups. For all layups, except the [0] 8 , the yielding of the metal 
matrix caused the fiber stress concentration factor to increase with 
increasing load. This is contrary to the behavior seen in 
homogeneous materials where yielding causes the stress concentration 
to drop. For the [0] 8 laminate, yielding of the matrix caused a 
decrease in the fiber stress concentration. The finite element 
analysis predicted these trends correctly. 
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I. INTRODUCTION 


Metal matrix composites (MMC) are being considered for many 
applications in advanced aerospace structures primarily because of 
their superior dimensional stability and high temperature capability 
characteristics compared to resin matrix composites. In addition, 
continuous fiber-reinforced MMC offer very high stif fness-to-weight 
and strength-to-weight ratios when compared to monolithic metals. 

Most of the projected aerospace applications will be required to be 
damage tolerant. In order to perform a damage tolerance analysis, it 
is important that the stress state around notches be well defined. 

In continuous fiber-reinforced MMC, the matrix can yield while the 
fiber remains almost perfectly elastic. This combination of 
constituent properties can lead to stress concentrations at slit tips 
that are different than for elastic or elastic-plastic homogeneous 
materials. The implications of the matrix yielding have not been 
widely recognized or understood in the materials community. The 
purpose of this paper is to examine experimentally and analytically 
the matrix plasticity effect around slits, in particular, the effect 
on fiber stress concentrations. 

Johnson, Bigelow, and Bahei-El-Din [1], and Poe and Sova [2] 
conducted a number of fracture tests on boron/aluminum (B/Al) 
laminates containing slits. Fiber and matrix damage was monitored by 
using radiography. Johnson, Bigelow, and Bahei-El-Din [1] found that 
first damage could be predicted by analytically monitoring the stress 
state in the boron fibers adjacent to the slit. When the fiber 
stress reached a predetermined value, damage occurred. First fiber 
failure occurred at approximately one half of the specimen's ultimate 
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strength for the slit unidirectional specimens. The first fiber 
failure in specimens with slits for [±45] 2s , [0/±45] s , and [0 2 /±45] s 
laminate orientations occurred at or very near the specimens ultimate 
strength. Since the damage does not usually propagate in a self 
similar manner in aluminum MMC, traditional fracture mechanic 
quantities such as stress-intensity factors may be questionable for 
general application to laminated metal matrix composites. Failure 
predictions using the stress (or strain) state in the load carrying 
fibers just ahead of the slit may be the most viable approach for 
assessing damage tolerance [1-3]. 

This paper will present experimental data for silicon-carbide/ 
aluminum composites in a number of different layups. The strain 
state ahead of the slit tip was determined experimentally with small 
strain gages. The strain distribution ahead of the slit will be 
presented for load levels that result in material responses ranging 
from elastic to highly elastic-plastic. The stress-strain state will 
also be predicted analytically using a three-dimensional finite 
element program. The predictions will be compared to the 
experimental data to assess the suitability of the analysis model. 


II. EXPERIMENTS 

A. Specimens, Materials, and Preparations 

Continuous fiber silicon-carbide/aluminum composite laminates 
with silts were tested statically to failure. The composite matrix 
material was 6061 aluminum which was annealed prior to testing. The 
fibers were 0.14-mm diameter silicon-carbide with a surface coating 
suitable for aluminum matrix applications. The fiber, designated 
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SCS 2 , and the composite were manufactured by AVCO Specialty Materials 
Division, Lowell, MA, in the form of 450 by 450 mm plates.* Several 
of the supplied composite plates were slightly bowed. This resulted 
in some of the specimens being slightly curved. This problem will be 
addressed later in the text. The supplied laminates contained fibers 
that had lower moduli than suggested by the manufacturer as explained 
in reference [4]. The composite constituent properties are given in 
Table I while the laminate elastic moduli are given in Table II. 

Five different layups were examined: [0] 8 , [0 2 /±45] s , [0/90] 2s , 

[0/±45/90] s , and [±45] 2s . Each specimen was 1.6 mm thick and was cut 
into a 51 by 164 mm rectangle using a diamond wheel saw. Either a 
9.5 or 19 mm long center slit was electronically discharged machined 
(EDM) into each specimen. The slits were approximately 0.3 mm wide. 

B. Experimental Procedures 

A row of five small strain gages (each approximately 1 mm 
square) was bonded just ahead of the slit tip on one side of the 
specimen to measure axial strain. The strain gages were spaced 2.0 
mm apart, centerline to centerline, as shown in Figure 1. 

The specimens were loaded in a tensile displacement control mode 
using a servo-hydraulic test machine. Approximately 39 mm on each 
end of the specimen was used for gripping. A computer was used to 
record the strain versus load data for each strain gage. The 
specimens were loaded at approximately 0.5 mm per minute until they 
failed. 

* Trades names are used to adequately describe the materials used 
and their use does not imply an endorsement by the National 
Aeronautics and Space Administration. 
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C. Results and Discussions 

The failure stresses of the slit specimens are given in Table 
II. The strain gage results are shown in Figures 2 through 6. Each 
figure shows the data for both the 9.5 and the 19.0 mm slit lengths 
in a given laminate. The normalized strain distribution ahead of the 
slit tip is given for a low load level and a high load level. The 
strain distribution at the low load level is considered to be almost 
entirely elastic, while the higher stress level causes extensive 
yielding in the matrix material. The measured strains are normalized 
by the calculated remote elastic strain for each given load level. 

The calculated remote elastic strain at a given load was found by 
dividing the applied stress by the calculated longitudinal modulus of 
the laminate given in Table II. The actual remote strain may be 
nonlinear at the higher load levels and would, therefore, be greater 
than the calculated remote elastic strain. 

As mentioned earlier, some of the test specimens had a slight 
curvature in the longitudinal direction. Therefore, upon tensile 
loading, some bending stresses were introduced. In order to 
compensate for this bending effect, the strain distributions were 
shifted such that the average value of the normalized strain for the 
elastic stress level was 1.6 for the 19-mm slit and 1.23 for the 
9.5-mm slit. The values 1.6 and 1.23 are the ratios of the gross 
cross-sectional area to the net cross-sectional area for the 
respective slit lengths. For equilibrium, the same load must pass 
through both the net and gross sections. Thus, the average elastic 
strain in the net section will be higher than the remote strain, by 
the factor of the ratio of the gross area to the net area. 
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The strain concentrations just ahead of the slit tip (at 0.5 mm) 
are higher for the plastic strain distribution than for the elastic. 
Even if the strain data curves were to be shifted vertically such 
that they came together at 8.5 mm from the slit tip, the plastic 
strain concentration would still be higher than the elastic. The 
only exception is the [0]g laminate. Unidirectional laminates 
develop long, narrow plastic zones at the slit tip due to the intense 
shear stresses. These zones are parallel to the fibers as noted by 
Jones and Goree [5] and shown in Figure 7. This shear yielding 
greatly reduces the strain concentration ahead of the slit. 

Therefore, in Figure 6, the strain distribution for the [0]g is 
rather flat for the high load, plastic case and the slit-tip strain 
concentration is lower for the plastic case than for the elastic 
case. 

For all laminates containing 0° plies, the axial strain measured 
at a point on the surface should be a good indicator of the stress in 
the 0° fibers beneath. The experimental data clearly indicate that 
the yielding of the matrix material causes the strain concentration 
near the slit tip to rise for all layups examined except the 
unidirectional case. Therefore, for all lavuos except the 
unidirectional, the axial stresses in the 0^ fibers iust ahead of the 
slit are significantly higher due to the plastic deformation of the 
matrix material than would be predicted with an elastic analysis. 

The analytical prediction capabilities will be examined in the next 
section. 


6 




III. ANALYSIS 


A. Homogeneous Materials 

In homogeneous metallic materials such as aluminum alloys, 
yielding reduces the slit-tip stress concentration, as illustrated in 
Figure 8. Here the stress distribution ahead of the slit, calculated 
using a finite element program, is shown for a center cracked, 
aluminum specimen. The specimen is 19.26 mm wide, 84.07 mm long, 
with a slit of length 5.49 mm located in the center of the specimen. 
The aluminum material is modeled using the same properties as the 
matrix material in the composite laminates. Predictions based on a 
stress failure criteria, using the elastic response, would be 
conservative for a specimen that had experienced yielding. However, 
this is not, in general, the case for metal matrix composites. As 
was shown by the experimental results of the previous section, 
yielding causes the slit-tip strain concentration to increase. This 
will be confirmed analytically in the following section. 

B. Metal Matrix Composites 

In metal matrix composites, the matrix yields, whereas the 
fibers generally behave elastically until they break. Therefore, to 
model the behavior of metal matrix composites, the elastic-plastic 
behavior of the matrix and the elastic behavior of the fiber must be 
accounted for. 

1. Finite element analysis 

An analysis that models this two phase behavior was used in the 
present study. This analysis was conducted with a three-dimensional 
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finite element program [6] called PAFAC (Plastic And Failure Analysis 
of Composites) , which was developed from a program written by 
Bahei-El-Din et al. [7,8]. PAFAC uses a constant strain, 
eight-noded, hexahedral element. Each hexahedral element represents 
a unidirectional composite material whose fibers can be oriented in 
the appropriate direction in the structural (Cartesian) coordinate 
system. 

Material model . - The PAFAC program uses a continuum material 
model developed by Bahei-El-Din and Dvorak [8,9,10] to represent the 
essential aspects of the elastic-plastic behavior of composite 
laminates. The material model is briefly described in the Appendix. 

Description of meshes . - Previous analyses using PAFAC [1,6] 
used a radial mesh, such as shown in Figure 9(a). However, at the 
time of the previous analysis [1,6] the authors felt that a high 
shear zone did not exist for slit laminates containing both 0° and 
angle plies. Recently, Post et al. [11] found rather significant 
shear zones at a slit tip in a [0/+45] s layup of boron/aluminum. 
Figure 10 shows a typical Moire fringe pattern found by Post et al. 
which illustrates the high shear band at the slit tip. It was felt 
that a radial mesh would not adequately model these high shear zones 
typically seen in front of the slit tip in 0° layers. Thus, a 
rectangular mesh, with a rectangular slit tip, as shown in Figure 
9 (b) , was also used in the present work to determine the effects of 
mesh type. Both meshes shown in Figure 9 have the same size element 
at the crack tip. The crack tip element was sized to represent one 
fiber spacing (0.178 mm). 
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2. Comparison of experimental and analytical results 

Figure 11 compares the predictions using the rectangular and 
radial meshes with the experimental results for the [0] 8 laminate 
with 2a = 19 mm. As before, these results are presented as 
normalized strains ahead of the slit tip. The solid and dashed 
curves show the predictions using the rectangular mesh for remote 
applied stresses of 12 and 400 MPa, respectively, while the two 
dash-dot curves represent results from the radial mesh at the same 
stress levels. Experimental results are shown by the symbols. In 
all cases, the lower stress level represents an elastic case and the 
higher stress level represents an elastic-plastic case. For this 
layup, the yielding of the matrix causes the slit-tip strain 
concentration to decrease rather than increase (as shown in Figure 
6) and the analysis using the rectangular mesh predicts this trend 
correctly. The analytical predictions using the radial mesh 
incorrectly predicted that the strain concentration would increase 
with increasing plasticity. For both the elastic and elastic-plastic 
stress levels, the rectangular mesh produced predictions that were in 
better agreement with the experimental results in the region close to 
the slit tip. The elastic-plastic prediction crosses over the 
elastic prediction very close to the experimental data. Thus, the 
rectangular mesh will be used for the remainder of the analyses 
presented. 

Figure 12 compares the analytical predictions using the 
rectangular mesh only to the test results for the [0] 8 laminate with 
2a = 19 mm. The solid curve and the circular symbols represent the 
normalized strains for an applied stress of 12 MPa and the dashed 
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curve and the square symbols are for an applied stress of 400 MPa. 

For this layup, the analytical predictions and the experimental 
results are in very good agreement for both the elastic and 
elastic-plastic cases over the entire region. 

Figure 13 compares the analytical predictions with the 
experimental results for the [0/+45/90] s layup with 2a = 19 mm. The 
two curves represent the analytical predictions while the symbols 
represent the experimental data. The results are shown for applied 
stress levels of 7 and 100 MPa. There is good agreement between the 
experimental results and the analytic predictions. Here the analysis 
shows that the slit-tip strain concentration factor increases with 
increasing plasticity, as was shown experimentally for this layup. 

Figure 14 compares the analytical predictions with the 
experimental results for the [0 2 /±45] s with 2a = 19 mm. The two 
curves represent the analytical predictions while the symbols 
represent the experimental data. The results are shown for applied 
stress levels of 11 and 200 MPa. Again, we see good agreement 
between the experimental results and the analytical predictions. 

The yielding of the matrix causes the strain concentration factor at 
the slit tip to increase for this layup. 

The strain distributions in the [0/90] 2s and the [±45] 2s layups 
were also predicted but are not shown here. The [0/90] 2s predictions 
were very close to the data. The analytical predictions did not 
agree well with the data for the [±45] 2s layup. In this case the 
strain gages indicated higher strain levels than predicted by the 
analysis. However, since the load carrying plies in this laminate 
are at 45°, the strain gages could not measure the strains in those 
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fibers. The finite element analysis could predict the fibers axial 
stresses and strains in the 45° PiY and found that the normalized 
strain concentration (and stress concentration) did increase with 
increasing plasticity. Further, the [±45] 2s layup may experience 
fiber rotation or extensive interlaminar yielding that the PAFAC 
program does not currently model, which would also result in low 
predictions from the analysis. 

Figure 15 shows the predicted longitudinal stress in the fiber 
closest to the slit tip, normalized by the applied stress, for all 
the layups containing 0° plies. The fiber stress concentration 
factor is shown for both an elastic and an elastic-plastic case for 
each layup. As previously noted, for the [0] 8 laminate, the plastic 
prediction is lower than the elastic prediction, while for all the 
other layups, the plastic predictions are higher than the elastic 
predictions. 


IV. DISCUSSION OF RESULTS 

We have shown experimentally and analytically that matrix 
yielding can significantly affect the fiber stress concentration. 
Basically two mechanisms of matrix yielding are working either to 
decrease the fiber stress concentration or increase it. 

The first mechanism was observed in the unidirectional layup. 
High shear stresses in the matrix material between the slit tip and 
the first continuous fiber form very long yield bands (as shown in 
Figure 7) that tend to isolate the remaining net section from the 
slit thus reducing the strain concentration due to the slit. This 
effect, to a lesser extent, may be present in any layup containing 0° 
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plies as shown by Post et al. (see Figure 10). 

The second mechanism is present in all MMC that experience 
matrix yielding; when the matrix yields the load is transferred to 
the load carrying 0° plies. The actual sequence of events is as 
follows: as the off-axis plies yield (and they yield before the 0° 

plies) , the load is transferred to the 0° plies; as the matrix in the 
0° plies yields, more load goes into the 0° fibers. Since for 
silicon-carbide/aluminum composites the matrix has approximately 20 
percent as much stiffness as the fibers, when the matrix loses 
stiffness, considerable load is transferred to the 0° fibers. The 
matrix is, of course, yielded more extensively closer to the slit 
tip. As a result, more load is transferred locally into the fibers 
in that region, thus, increasing the strain concentration at the slit 
tip. 

V. CONCLUSIONS 

Continuous fiber silicon-carbide/aluminum composite laminates 
with slits were tested statically to failure. Five different layups 
were examined: [0] 8 , [0 2 /±45] s , [0/90] 2s , [0/±45/90] s , and [±45] 2s . 

Either a 9.5 or a 19 mm slit was machined in the center of each 
specimen. The strain distribution ahead of the slit tip was found 
experimentally with a series of strain gages bonded ahead of the slit 
tip. A three-dimensional finite element program (PAFAC) was used to 
predict the strain distribution ahead of the slit tip for several 
layups . 

For all layups, except the [0] 8 , the yielding of the metal 
matrix caused the fiber stress concentration factor to increase with 
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increasing load. This increase is primarily due to the load being 
transferred out of the yielded matrix material and into the load 
carrying 0° fibers. The finite element analysis predicted this trend 
correctly. This is contrary to the behavior seen in homogeneous 
materials where yielding causes the stress concentration to drop. 

For the [0] 8 laminate, yielding of the matrix caused a decrease 
in the fiber stress concentration. This decrease was due to 
shear stress bands that form at the tip of the slit and grow 
parallel to the fibers. These shear band essentially isolate the 
notch effect from the remaining net section, thus lowering the strain 
concentration. Again, the analysis agreed with the experimental 
results for this layup. 
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APPENDIX 

MATERIAL MODEL OF PAFAC PROGRAM 

The model consists of an elastic-plastic matrix unidirectionally 
reinforced by continuous elastic fibers. Both constituents are assumed to be 
homogeneous and isotxcpic. The fibers cure assumed to have a very small 
diameter, so that although the fibers occupy a finite volume fraction of the 
composite, they do not interfere with matrix deformation in the two transverse 
directions, but only in the axial (fiber) direction. Figure 16 shows a 
schematic of this lamina model. It can also be represented by parallel fiber 
and matrix bars or plates with axial coupling, as illustrated in Figure 16. In 
Figure 16, the fiber (axial) direction is parallel to the x 3 -axis, and the x^- 
and x 2 ~axes represent the transverse directions. 

If the Cartesian coordinates are chosen so the x^ coincides with the fiber 
direction, the second-order tensors of the independent stress and strain 
components, a and e, are expressed as 

a = [a ll ° 22 a 33 ° 12 ° 13 ° 
e = 11 £ 22 £ 33 7 12 1 13 ^ 2 b 

where -y „ = 2e ^ (i, j = 1, 2, 3; i * j) are the engineering shear strain 
components. 

For equilibrium and compatibility, several requirements are imposed on the 
material model shewn in Figure 16. The stress average in each constituent can 
be related to the overall composite stress a in the axial (fiber) direction as 
follows: 
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a 33 V f (a 33 J f + V m (<7 33^m 


A bar over a symbol indicates overall composite stress or strain, and the 

subscripts f and m denote quantities related to the fiber and matrix. The 

volume fractions v~ and v are such that v- + v =1. The other stress 

r m r m 

components in each constituent were assumed to be uniform and to obey the 
following equilibrium equations: 

a ll = ^ll^f 52 (a ll*m 
*22 = <*22>f" ^22>m 
a 12 = ^12* f = (a 12^m 
a 13 = * a l3*f “ (<7 13? m 
a 23 = (ff 23* f = (<T 23* m 

The only constraint in the model is in the axial (fiber) direction; the 
matrix and fiber must deform equally. Thus, 

e 33 = 3^ f = (e 3^ m 

The other strain components can be related to the overall strain, 7 as 
follows: 

r ij = v f (£ i^ f +v !l £ ijm <ij*33> 

Since the fibers are elastic up to failure, the inelastic strains of the 

lamina are caused by matrix deformation. Because the fiber imposes an elastic 
constraint on the matrix which affects the shape of the lamina yield surface, 
additional kinematic components appear in the hardening rule of the lamina and 
influence the magnitude of the overall plastic strains. All aspects of the 
yield behavior were examined and accounted for in the formulation of the lamina 
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constitutive equations. These equations are explicitly described in [6,8]. The 
stress-strain curve of the matrix material was modeled with a Ramberg-Osgood 
equation [10] . 


16 



REFERENCES 


[1] Johnson, W. S.; Bigelow, C. A.; and Bahei-El-Din, Y. A.: 
Experimental and Analytical Investigation of the Fracture 
Processes of Boron/Aluminum Laminates Containing Notches. NASA 
TP-2187, National Aeronautics and Space Administration, 
Washington DC, 1983 

[2] Poe, C. C. , Jr. and Sova, J. A.: Fracture Toughness of Boron/ 
Aluminum Laminates With Various Proportions of 0° and +45° 
Plies," NASA TP-1707, National Aeronautics and Space 
Administration, Washington DC, November 1980. 

[3] Poe, C. C. , Jr.: A Unifying Strain Criterion For Fracture of 
Fibrous Composite Laminates, Engineering Fracture Mechanics . Vol 
17, No. 2, pp. 153-171, 1983. 

[4] Johnson, W.S.; and Wallis, R.R.: "Fatigue Behavior of 
Continuous-Fiber Silicon Carbide/Aluminum Composites," Composite 
Materials: Fatigue and Fracture . ASTM STP 907, H.T. Hahn, Ed., 
American Society for Testing and Materials, Philadelphia, 1986, 
pp. 161-175. 

[5] Goree, J.G.; and Jones, W.F. : Fracture Behavior of 
Unidirectional Boron/Aluminum Composite Laminates. NASA 
CR-3753, National Aeronautics and Space Administration, 
Washington DC, December, 1983. 

[6] Bigelow, C. A.; and Bahei-El-Din, Y. A.: Plastic and Failure 

Analysis of Composites (PAFAC) . LAR-13183, COSMIC, Univ. of 

Georgia, 1983. 

[7] Bahei-El-Din, Y. A.; Dvorak, G. J. ; and Utku, S.: Finite 

Element Analysis of Elastic-Plastic Fibrous Composite 
Structures. Comput. & Struct., vol. 13, no. 1-3, June 1981, pp. 
321-330. 

[8] Bahei-El-Din, Y. A.: Plastic Analysis of Metal-Matrix Composite 

Laminates. Ph.D. Disser. , Duke Univ., 1979. 

[9] Dvorak, G. J.; and Bahei-El-Din, Y. A.: Plasticity of Composite 

Laminates. Research Workshop on Mechanics of Composite 
Materials, Duke Univ., Oct. 1978, pp. 32-54. 

[10] Bahei-El-Din, Y. A.; and Dvorak, G. J. : Plastic Yielding at a 

Circular Hole in a Laminated FP-Al Plate. Modern Developments 
in Composites Materials and Structures, J. R. Vinson, ed., 
American Soc. Mech. Eng., c. 1979, pp. 123-147. 


17 


[11] Post, D. ; Czarnek, R. ; Joh, D. ; Jo, J.; and Guo, Y: 
Elastic-Plastic Deformation of a Metal-Matrix Composite Coupon 
with a Center Slot. NASA CR-178013, National 

Aeronautics and Space Administration, Washington DC, November, 
1985. 

[12] Ramberg, W. ; and Osgood, W. R. : Description of Stress-Strain 

Curves by Three Parameters. NACA TN 902, 1943. 


18 


/ 


Table I. Composite Constituent Mechanical Properties 


Elastic modulus, GPa 
Poisson* ratio 


SCS 2 fiber 
0 . 14 mm diameter 

340 

0.25 


6061-T0 Aluminum 
72.5 
0.33 


Table II. Failure Stresses and Calculated Elastic Moduli 


Laminate 

Slit length, mm 

Failure Stress, MPa 

E L , GPa 

[0] 8 

9.5 

793 

190 


19 

545 


[0 2 /±45] s 

9.5 

363 

154 


19 

258 


[ 0/90 ] 2s 

9.5 

269 

153 


19 

200 


[ 0/ ±4 5/90 ] s 

9.5 

259 

137 


19 

200 


[±45] 2s 

9.5 

177 

118 


19 

88 
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Figure 1. Specimen configuration. 
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Distance from the slit tip, mm 

Figure 2. Normalized strain in front of slit tip for [0/90] 2s layup. 
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Figure 3. Normalized strain in front of slit tip for [02/±45] s layup. 
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Figure 4. Normalized strain in front of slit tip for [0/±45/90] s layup. 





25 


Distance from the slit tip, mm 

Figure 6. Normalized strain in front of slit tip for [0] s layup. 
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Figure 7. Brittle lacquer crack pattern at 12,500 N[3]. 




Distance from crack tip, mm 

Figure 8. Stress distribution ahead of slit tip in homogeneous aluminum. 
2a = 5.86 mm. 
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Figure 9. Finite element meshes 
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Figure 11. Comparison of rectangular and radial mesh results for [0 ]q 
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Distance from slit tip, mm 



p with rectangular mesh 
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Figure 15. Fiber stress concentrations for layups containing 0°plies. 
2a = 19 mm. 
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